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Abstract 

The bound state solution of the radial Schrodinger equation with the generalized Woods-Saxon 
potential is carefully examined by using the Pekeris approximation for arbitrary (. states. The 
energy eigenvalues and the corresponding eigenfunctions are analytically obtained for different n 
and t quantum numbers. The obtained closed forms are applied to calculate the single particle 
energy levels of neutron orbiting around ®®Fe nucleus in order to check consistency between the 
analytical and Gamow code results. The analytical results are in good agreement with the results 
obtained by Gamow code for £ = 0. 
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I. INTRODUCTION 


The Woods-Saxon potential was firstly proposed by R.D. Woods and D.S. Saxon in 
order to explain the elastic scattering of 20 MeV protons by medium and heavy nuclei 


approximately sixty years ago[l|. Thereafter the Woods-Saxon potential has taken a great 
deal of interest over the years and has been one of the most useful model to determine 

nn riQ 

the single particle energy levels of nuclei[2H^ and the nucleus-nucleus interactions [5H7|. 
The modified version of the Woods-Saxon potential consists of the Woods-Saxon and its 


derivative called the Woods-Saxon surface potential and is given by 

V (r) = 


id, 


( 1 ) 


where Vq and Wq represent the depths of the potential well. R and a are the radius of 
the potential and the width of the surface diffuseness, respectively. In Fig{T], a form of the 
generalized Woods-Saxon potential (GWS) versus the internuclear distance is shown for the 
given potential parameters with several Wq values. The surface term in the generalized 
Woods-Saxon potential induces an extra potential pocket, especially, at the surface region 


of the potential, and this poc 
of some nuclear reactions 


^t is very important in order to explain the elastic scattering 
Moreover, the Woods-Saxon surface potential induces a 
potential barrier for Wq < 0 so that it could be used in explaining of the resonant states 
(quasi-bound states) in nuclei. There are some special cases of the generalized Woods-Saxon 
potential: The GWS potential is reduced to the standard Woods-Saxon form for Wq = 0 


and the square well potential for Wq = 0 and a —)■ 0 n|. Furthermore the GWS potential is 


reduced to the Rosen-Morse potential 


12| for R = 0 


m- 


The relativistic treatment of a Dirac particle in the Woods-Saxon potential field is ex¬ 


amined in three dimensions for £ = 0 


14l | . Moreover, the transmission coefficient and bound 




state solutions of one dimensional Woods-Saxon potential are analytically studied[l5j. The 
Klein-Gordon equation in the presence of a spatially one-dimensional Woods-Saxon poten¬ 
tial is also examined. The scattering state solutions are obtained in terms of hypergeomet¬ 
ric functions and the condition for the existence of transmission resonances is derived jl^. 
Furthermore s-wave solution of the Dirac equation for a particle moving in the spherically 


symmetric Woods-Saxon poten 


symmetry limit is examined 


ial under the conditions o 


17| | and is discussed in Refs. 


the exact spin and the pseudospin 


18 


19|. 
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It is known that the exact analytical solutions of the wave equations (Schrodinger, Dirac, 
etc.) are very important since a closed form of the wave function is more convenient than 
the wave function obtained by numerical calculation in explaining the behavior of the sys¬ 
tem under consideration. Unfortunately, there are few potentials such as harmonic oscilla¬ 
tor, Coulomb and Kratzer potentials efc.[n| which have the exact analytical solution with 
centrifugal term. In literature, there are some effort about obtaining the approximate an¬ 
alytic solutions of the wave equations in terms of the £ 7 ^ 0 case: The most widely used 
approximation is introduced by Pekeris[^ for the exponential-type potential so that this 
approximation is based on the expansion of the centrifugal barrier in a series of exponen¬ 
tials depending on the internuclear distance up to second order. Greene and Aldrich 21 ] 
proposed another approximation for the centrifugal term 1/r^ However, 


this approximation is valid only for small values of the screening parameter 6 


22 |. 


Recently, the GWS potential has been examined by using the Nikiforov-Uvarov (NU) 


method 


23| for i = 0 state However, in the paper authors have obtained the energy 


eigenvalue equation as R independent due to r — R = r transformation so that the poten¬ 
tial is reduced to the Rosen-Morse potential 13(|. It has been noted that the analytical and 
the numerical results are inconsistent for the GWS potential with i = 0 state 13]. Simi- 

n nn 

lar results with Ref.llOH can be found in Refs. 24-27] for the relativistic or non-relativistic 

n 

wave equations. In Ref. [241] . the approximate analytical solution of the Schrodinger equa¬ 
tion for the standard Woods-Saxon potential is obtained for any i state by using the NU 


method. In Ref. 


equations. In Ref. [26 


25Lt 


le GWS potential is examined for the Klein-Gordon and Schrodinger 


271 ] the Woods-Saxon potential has been analyzed for both the radial 


Schrodinger and Klein-Gordon equations by using the Pekeris approximation. The authors 


in Refs. 


26 


27j ] have used z{r )=—transformation and have obtained R dependent eigen- 


l+e 


value equation by using the Nikiforov-Uvarov method [26l. l27l ]. However, we have realized 
that the Nikiforov-Uvarov method can not take into account the boundary condition cor¬ 
rectly since the Woods-Saxon potential has different character close to r=R. Therefore, in 
this article, we have carefully examined the radial Schrodinger equation with the generalized 
Woods-Saxon potential by using the Pekeris approximation in terms of the correct boundary 
conditions. In the next section, we present the calculation procedure. Then, in section IHII 
is devoted to the summary and conclusion. 
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II. THE ENERGY EIGENVALUES AND EIGENFUNCTIONS 


The generalized Woods-Saxon potential or special forms of it are very useful in order 
to describe the interactions between the systems, especially, in nuclear physics. In order 
to explain the single particle energy levels or elastic scattering of nuclei, the Woods-Saxon 
potential is generally used. Since the interactions usually occur at the surface region of nuclei 
for both bound and continuum states, the form of the potential at the surface is crucially 
important. Therefore the surface term in Eq.([T]) would be a very convenient model in order 
to calculate the single particle energy levels of nuclei. When we consider two-body system 
with the reduced mass fi moving under the generalized Woods-Saxon potential, the effective 
potential is. 


VefAr) = U(r) + W(r) 


l + e^~ + ^ 2/ir2 


( 2 ) 


where u = . rrin and ruA are the atomic mass of the neutron and the core nucleus re- 

rrin+mA 

spectively. There is no analytical solution of Eq.(|2]) due to polynomial form of the centrifugal 


barrier term. Therefore, we have to use an approximation for the centr ifug al 


;erm similar to 


other authors 26|]. In literature there are few approximation procedure 20|, l2]|. One of them 


is the Pekeris approximation 2^ based on an expansion of the centrifugal barrier depending 


on the internuclear separation up to second order 


26]. 


The quasi-analytical solution of the effective potential in Eq.(l2]) with the Pekeris ap¬ 
proximation can be obtained within the framework of the Nikiforov-Uvarov (NU) or 
asymptotic iteration (AIM) methods as follows. 


nr(nr + l) —/3^—J2 + (l + 2nr)e + 2e^ + (l + 2nr + 2e)-\^e‘^ + — /3'^ = 0, = 0,1, 2.... (3) 


with the following dehnitions: 

2fia^(E - 6Co) 


e^ = 


2 2/ua^(Vo - SCi] 








2 2iiar5C2 
7i = 




,1 = (4) 


a 






C 2 = 


a 




48 




R 


a = 


5 = 


In Ref. 13| , the consistency of the analytical results of Ref. lOj] with the Gamow code 3l|] has 


h? 

2 pR2 


been checked by calculating the single particle energy levels of the neutron moving around 
the ®®Fe nucleus for the given potential parameters. As can be seen in Ref, jl^ . the results 


are inconsistent with the numerical calculations. We have also conhrmed that there are 
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inconsistencies between Egand the Gamow code results for £ = 0. In Eq.([3]), if one uses 


Wo = 0, the results of Ref. 


26| can be obtained for arbitrary i states. Furthermore if we 


take £ = 0 in Eq.([n]), we get the results of Ref. jl0|. Therefore we should say that Eq.dH]) 
determined by the NU method is incorrect. If we used any analytical solution methods 
such as the asymptotic iteration method (AIM)[28|, the supersymmetry (SUSY)[l2|, etc. to 
solve the corresponding equations with the generalized Woods-Saxon potential, we would 
End same results in Eq.(|3]). In literature, there are similar calculations for the analytical 
solution of the generalized Woods-Saxon potential with the relativistic or non-relativistic 


wave equations by using the analytical methods 


29l |. The origin of the problem is due to the 


boundary conditions so that the analytical methods cannot take into account them correctly 
since the Woods-Saxon potential has different characteristic neighborhood r=R. Therefore 
we carefully have to examine the boundary conditions. In order to get the asymptotic 
behavior of the wave function Un£{z), we can use the Nikiforov-Uvarov method 23|] and easily 
get (j){z) = z^{l — . As a result, the asymptotic solution of the wave function is. 


Untiz) = Z%1 - zffntiz), 


(5) 


where z=—and The wave function in Eq.([5]) satishes the boundary 

1+e a 

conditions, he., Uni{r ^ 0,;^ ^ 1) —)■ 0 and Une{r oo, z ^ 0) ^ 0. The Schrodinger 
equation becomes for Eq.([5]), 

+[l+2e — {2 + 2e+2r])z] ^ ^ ^ “[“( 7 i+ 72 )+^+h + (^+h)^]/n£(^) =0- (6) 


dz"^ 


dz 


It is known that the hypergeometric equation 30(] is defined as 


.d^wiz) . , , .dwiz) , / s 

z{l — z) — j ^ ^ + [c — {a + b + l)z]—— - abw{z) = 0, 


dz'^ 


dz 


(7) 


and one of the solutions is w{z) = 2 Fi{a,b;c; z)^^. In order to get a,b,c parameters we 
compare Eq.([S]) with Eq.([7]) and find. 


a = -(1 T ^ 1 + 47? + 472 + 2 e + 2 r/), 
6 = ^(1 ± + 472 + 472 + 2 e + 2r]), 


c = 1 + 2e. 


( 8 ) 


Consequently we have Une{z) = Y(1 — z)^ 2 Fi{a,b-, c; z). To study in the vicinity of 2 ; = 1, 
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Q. 


we use the relation 


1^/1 \ r(c)r(c —a —6) , , 

2 Fi(a, 6;c;|/) = —-—--+ h- c+l]l-y) 


r(c-a)r(c-6)- 
r(c)r(a + 6-c) 
r(a)r(6) 


(1 - yY “ ^ 2 -Fi(c - a,c-h]c- a - h + l]l- y). (9) 


2 -Fi(a, 6; c; 0) is equal l[30|. Therefore, by using Eq.([9]) and the boundary condition 
0, ^ ^ 1) = 0, we get 

r(a + h — c) r(c — a) r(c — b) 




( 10 ) 


r(c-a-6) T{b) r(a) 

where rj = iX and A = \/(3“^ — 'Ji — Evaluating Eq. fllOp for the given a,b,c parameters in 
Eq-dH]), we obtain 

r( 2 iA) r(A + 1^1 + 4^2 ^ g _ r(l - 1^1 + i'yf + 47 I + £ - if]) 

r(-2iA) r(A - 1^1 + 47 ^ + 47 I + e + ir]) r(| + + e + irii) 

x(i + gi?/a) 2 ^ _ _1_ 

In Eq. flllj) . S> 1 for a given realistic parameters. Therefore we can use an approximation 
1 _j_gH/a ^ gij/a small errors in the eigenvalues. We can easily get the following equation 
by using relation, 


2i [argr{2iA)-argr{i + i^l+472+47|+e+i?7)-argr(i-i-^l+472+47|+£+j)?)+.^] _ 


( 12 ) 


Therefore, the corrected energy eigenvalue equation in a closed form becomes, 

argT{2iX)-argT{^ + ]^^Jl^^Ail^^^ + e + ir|) (13) 

-argV{]- - I + A-fl + + e + ig) + — = {ur + ^)7r, Ur = 0, ±1, ±2,..., 


where is the radial node number. The quantum number is n=nr + 1. In order to test 
the accuracy of Eq. flT^ . we calculate the single particle energy levels of neutron rotating 
around ^®Fe nucleus by using the potential parameters given in Ref, [l^ . The Woods-Saxon 
potential parameters are Vo=40.5 + 0.1374=47.78 MeV, R=4.9162 fm and a=0.6 fm. Here 
A is the atomic mass number of ^®Fe nucleus. The reduced mass consists of neutron mass 
rrin = 1.00866m, and ®®Fe core mass = 56m. In Table [H we show agreement between 
the energy eigenvalue equation given by Eq. (ll3l) and the numerical calculation obtained by 


Gamow code 


3l| for the neutron plus ^®Fe nucleus system with several rir quantum numbers 


and ITo parameters. There are small inaccuracy in the analytic and numeric results since 
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we have made the approximation in Eq. ifTT]) . It might be seen that neutron is unbound for 
fir = 3, Wo = 0 and Eq. ffT^ gives En^e = 62.9775MeV, but this Ending is not acceptable 
as physically and is only a mathematical result which satisfies Eq. flldp . There are similar 
situations for = 3, Wq = 50MeV; Hr = 3, Wq = lOOMeV; and rij. = 2, 3, Wq = —SOMeV. 
However there are the quasi-bound states (resonant states) for = 2,Wo = —100 MeV 
since the nuclear potential has a potential barrier inducing the resonant states in Fig{T] Our 
result for = 2, Wq = —100 MeV is physically incorrect. In order to calculate the en ergy 
eigenvalues of the resonant states, the Complex Scaling Method(CSM) can be used|32j. 
Another interesting point is that we cannot calculate the bound state energy eigenvalues by 
using Eq. flT^ for riy. = 2,3,Ho = lOOMeV since the right and left sides of Eq. dT^ do not 
intersect in real line. It should be noted that the Astate solution of the generalized Woods- 
Saxon potential in terms of the Pekeris approximation is valid only for small a values. 

The radial wave function corresponding to the eigenvalue equation Eq.([T3]) in terms of 
Eq.(j4]) and Eq.(jH]) can be written in a closed form as follows, 


Unt{r) = N 


l + e~ 


1 - 


1 + e" 


(14) 


X 2Fi{- - 


- ^^1 + 472 + 472 + £ + r;, 1 + + 472+ 472 + e + r]]l + 2 e ] -—^^), 

1 “t” e ^ 


'2 2 

where N is the normalization constant. The unnormalized wave function fulfilling the bound¬ 
ary conditions is shown in Figj2]for several radial quantum numbers n. 


III. CONCLUSION 

We have studied the approximate analytical solution of the Schrodinger equation in the 
presence of the generalized Woods-Saxon potential by using the Pekeris approximation and 
the properties of the Hypergeometric functions for any £ states. We have seen that the 
Nikiforov-Uvarov method cannot take into account the correct boundary conditions for the 
generalized Woods-Saxon potential. Therefore we have carefully examined the asymptotic 
behavior of the wave function of the generalized Woods-Saxon potential and have obtained 
the corrected eigenvalue equation and the corresponding eigenfunction in the closed form 
for any £ states. We have also calculated the single particle energy level of neutron rotating 
around ^®Fe nucleus for the given potential parameters in order to check the consistency 
between the analytical and Gamow code results. We have shown that the obtained analytical 
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results in this study are in good agreement with the results obtained by the Gamow code 
for £=0 state. The resonant state solutions of the generalized Woods-Saxon potential in a 
closed form are in progress. 
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Wo{MeV) 

Ur 

^Analytical 


0 

0 

-38.3004 

-38.3002 

0 

1 

-18.2254 

-18.2227 

0 

2 

-0.2678 

-0.2663 

0 

3 

62.9775 

unbound 

50 

0 

-41.1965 

-41.1964 

50 

1 

-23.8789 

-23.8788 

50 

2 

-3.6472 

-3.6471 

50 

3 

52.0232 

unbound 

100 

0 

-45.4453 

-45.4446 

100 

1 

-29.1659 

-29.1642 

100 

2 

undetermined 

-7.8143 

100 

3 

undetermined 

unbound 

-50 

0 

-36.2136 

-36.2168 

-50 

1 

-12.8469 

-12.8504 

-50 

2 

18.5701 

unbound 

-50 

3 

64.2816 

unbound 

-100 

0 

-34.5902 

-34.5956 

-100 

1 

-8.0843 

-8.0902 

-100 

2 

19.2098 

17.74-hi(-8.36) 


TABLE I. Comparison of the analytical and numerical results for the single particle energy levels 
of neutron orbiting around ^®Fe nucleus with several Wq potential depth and rir quantum numbers 
for i = 0. The potential parameters are Vq = 40.5 + 0.13A = 47.78 MeV, R = 4.9162 fm and 
a = 0.6 fm. The reduced mass consists of neutron mass = 1.00866u, and ^®Fe core mass which 
is niA = 56u and its value is /j, = 0.990814u. 
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FIG. 1. Variation of the generalized Woods-Saxon potential as a function of the internuclear 
distance and several Wq values for Vo=40.5+0.13A MeV, R=4.9162 fm and a=0.6 fm. 
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FIG. 2. The unnormalized wave function of the generalized Woods-Saxon potential for 
Vb=40.5+0.13A MeV, 11^0=50 MeV, R=4.9162 fm and a=0.6 fm potential parameters with several 
n quantum numbers. 
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